Magnetic field dependence of pairing interaction in ferromagnetic superconductors 

with triplet pairing 
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It is developed a microscopic description of superconductivity in ferromagnetic materials with 
triplet pairing triggered by the exchange of magnetic fluctuations. Instead widely used paramagnon 
model we work with phenomenological spectrum of fluctuations in the orthorhombic ferromagnet 
with strong magnetic anisotropy. Depending of the field orientation parallel or perpendicular to the 
direction of spontaneous magnetization the effective amplitude of pairing interaction proves to be 
decreasing or increasing function of magnetic field that allows to explain the drastic difference in 
magnitudes of upper critical field in these directions. 

PACS numbers: 74.20.De, 74.20.Rp, 74.25. Dw 



I. INTRODUCTION 

Superconductivity in ferromagnetic state found in ura- 
nium compounds URhGe and UCoGe^ exhibits peculiar- 
superconducting properties such as extremely high up- 
per critical field" and the reentrant superconductivity 
in large external magnetic fields^,. These ferromagnets 
have orthorhombic structure with magnetic moment ori- 
ented along c axis. The observed upper critical field in 
the directions perpendicular to the magnetization proves 
much higher than the paramagnetic limiting field. This 
unexpected for the superconductors with triplet pairing 
behavior— is explained by the band splitting caused by 
exchange interaction 11 . On the other hand, the param- 
agnetic limitation in a superconductor with equal spin 
pairing for the field parallel to the spin quantization di- 
rection is absent^. However, experimentally, the upper 
critical field in the direction of spontaneous magnetiza- 
tion approximately coincides with the paramagnetic lim- 
iting field in these materials. The question, why the up- 
per critical field parallel to c axis is so much smaller than 
the H C 2 in the perpendicular directions, remains open. 

Another peculiar feature observed in UCoGe^— for the 
field directed perpendicular to the spontaneous magne- 
tization is a pronounced upward curvature. This fact 
cannot be explained by the field dependent increase of 
the effective mass probably responsible for the S-shape 
upper critical field along b axis behavior in the field re- 
gion from 5 to 15 Tesla. Indeed, the upward curvature 
of H C 2 (T) along a and b directions is observed already at 
much smaller fields where effective mass falls down with 
field increase 5 stimulating the opposite that is downward 
curvature tendency in the H C 2(T) behavior. The possi- 
ble explanation for anomalous upper critical field cur- 
vature pointed out in the paper— is attributed to the 
crossover between two phases in a two-band ferromag- 
netic superconductor s] 12 ' 13 . 

Here we discuss the circle of problems related to the up- 
per critical field behavior specific for superconducting fer- 
romagnets URhGe and UCoGe. Our treatment is based 
on semi-microscopic description of superconductivity in 



ferromagnetic materials with the Cooper pairing caused 
by the magnetic fluctuations. Instead of wide spread 
paramagnon or Fermi liquid approach to the nonphonon 
mechanism of superconductivity (see for instance^ and 
references therein) we shall use the phenomenological de- 
scription of magnetic fluctuations in the orthorhombic 
ferromagnet with strong spin-orbital coupling. Certainly 
the complete theory of triplet superconductivity in fer- 
romagnets has to take into account the frequency depen- 
dence of the pairing interaction. The invalidity of the 
Migdal theorem adds, however, the supplementary trou- 
bles to the Eliashberg type theory of superconductivity 
in heavy fermionic materials. Leaving this problem for 
the future investigations^ we shall see to what kind of 
qualitative conclusions one can come in frame of semi- 
microscopic weak coupling theory. 

To investigate the effect of the pairing interaction field 
dependence we drop out all the orbital effects. In other 
words we will solve the problem of the critical tempera- 
ture field dependence as if the magnetic field acts only on 
the electron spins. This simplification allows to demon- 
strate the pure effect of suppression of superconducting 
state for the field parallel to magnetization and the oppo- 
site effect of stimulation of superconductivity for the field 
directed in perpendicular to magnetization directions. 

It should be noted that the proposed mechanism of 
pairing interaction field dependence is not an alternative 
to the discussed in literature mechanism related to the ef- 
fective mass field dependence&i^— . Both of them make 
a contribution to the interplay of magnetism and super- 
conductivity in ferromagnetic compounds. The unified 
treatment taking into account both mechanisms in frame 
of some general approach is the problem for the future. 

The paper is organized as follows. The critical temper- 
ature field dependence for field orientations parallel and 
perpendicular to the direction of spontaneous magneti- 
zation are derived in two separate sections of the paper 
following each other. Both of them are based on the 
pairing interaction expressed through the components of 
susceptibility tensor derived in the Appendix. 
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II. PAIRING HAMILTONIAN 

Following the paper— we introduce the interaction be- 
tween electrons due to the spin fluctuations as 

Hint = ^ Di i ( k > k ')°"a5°il 7 a L a -k/9 a -k'7 a k'5- W 

kk'q 

The matrix is taken as the part of static spin suscep- 
tibility being an odd function of both of its arguments 

A,(~k,k') = A,(k, -k') = -A,(k,k'). 

It determines the pairing interaction for superconduct- 
ing state with triplet pairing in centrosymmetric crys- 
tal and has the structure corresponding to orthorhombic 
anisotropy in the spin space 

/ D x (k,k') \ 

Ai(k,k')= Aj(k,k') o . (2) 

V D z (k,k>) J 

The corresponding p-wave interaction matrix for 
isotropic Fermi liquid like liquid 3 He is 

D^(k,k') = -V(kk')6 l3 (3) 

The explicit form of matrices of interaction in the or- 
thorhombic ferromagnet are derived in the Appendix. 

After some algebra^ the hamiltonian Eqn.(l) can be 
rewritten as hamiltonian of pairing interaction 

Hint = q E ^/37i( k > k ') a L a -k£ a -k'7 a k'(5, (4) 
kk'q 

here 

V aj3l 8{k, k') = V ij (k,k'){i(Ti(jy) a s{i(jj<7y)^ i (5) 
I 



Vy(k,k') - ^TrDik^Sij - A,(k,k'). (6) 

The critical temperature is determined from the lin- 
earized self-consistency equation. In space homogenious 
case it has the following forrn^ 

A a/3 (k) = -T^^^ aAM (k,k') 

n k' 

xG X7 (k',u n )G^(-k',-io n )A jS (k'), (7) 

where G\ 7 (k' ,u> n ) is the normal metal Green function. 
The interaction functions -Dy (k, k') as well as the Green 
functions have the different form depending of the exter- 
nal field direction. We start from the simplest case of 
magnetic field parallel to the magnetization direction. 
III. MAGNETIC FIELD PARALLEL TO 
MAGNETIZATION AXIS 

For this field orientation the normal state matrix Green 
function is 

G = ( Gt M G^ = 1 
" \ G l J ' U iu n -^±^ B (h + H z y 

(8) 

where h is the exchange field and H z is the external field. 
Taking the matrix of order parameter as 

and substituting the matrices of the Green function and 
the order parameter into self-consistency equation after 
simple algebra we come to the system of equations 



At(k) = -T^^{ J D|(k,kOGtG^(k')+i?l(k,k')GtG^(k')}, 



(10) 



n k' 



A^(k) = -T^^{Di(k ) k')GlG|At(k') + £'|(k,k')GfG^(k')}, 



(11) 



n k' 



A^(k) = -It E { K (k ' k ' } ~ ^ (k ' k ' } 



Lr l Lr 2 "+" ( - t 1 Lt 2 



n k' 



A n (k')} • 



(12) 



Here we have introduced notations D±(k,k') = 
D x {k,k') ± D„(k,k') and G\ = &(\t?,u n ), C\ = 
G^(— k', — to n ) and similarly for the G\ and G\ Green 



functions. 

In the simplest case of superfluid 3 He with pairing in- 
teraction in the form ([3]) when D x — D y — D z we have 
three independent equations for AT, A^ and A T ^. More- 
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over, in absence of the exchange field h and the external 
field H z the Green functions for the Fermi particles with 
up and down spin coincide each other G^ — G"k The 
critical temperature has the same value for all the order 
parameter components A^, A^, A^, that corresponds 
to the phase transition to one of superfluid phases of 
3 He. The particular choice between them can be made 
by comparison of nonlinear terms in Ginzburg - Landau 
equation. 9 In presence of the external field the equality 
between the Green functions is violated ^ G"k It 
leads to the lifting of the critical temperatures degen- 
eracy for different components of order parameter. As 
result, the highest temperature phase transition occurs 
to the spin up-up pairing state known as A\ phase of 
superfluid 3 He. 

In our case of orthorhombic ferromagnet the equa- 
tions for order parameter components A^, A^, those are 
amplitudes of pairing of quasiparticles with spin up-up 
and spin down-down correspondingly, are coupled each 
other. The component A^ corresponding to pairing of 
quasiparticles from spin-up band with quasiparticles from 
spin-down band obeys the separate equation. In the fer- 
romagnetic state in view of the large band splitting the 
interband pairing superconducting state arises at defi- 
nitely lower temperature than the spin up-up and spin 
down-down states. In the paramagnetic state in general 
it is not true. Here we shall not be interested in the inter- 
band pairing superconductivity. So, we deal with system 
of two equations (flU)) and (TTTI) . 

All the interaction functions Di(k, k'), here i — x,y, z 
as one can see from Eqs. (IA15|1 . (|A16I) and (|A17|) have 
the same structure 



-^i(k, k ) — CLik x k x -\- b^kyky -\- Cjk z k z . 



(13) 



Hence, the system of equations ((TU)) and ([fi]) have three 
independent solutions with different critical tempera- 
tures 

At(k)=A, A+(k) =rr^; (14) 
At(k)=^ fcyj A+(k)=£^; (15) 

A*(k) = Ck z , A+(k) = C l h. (16) 

Let us assume that the largest critical temperature corre- 
sponds to solution given by Eqn. (|14[) . Then performing 
all necessary integrations and summation in Eqs. (|10[) and 
(fTTj) we come to the system of algebraic equations 



v f = (ajy + .gjV)A(T) 

^ = (ffj7? t + 9 J^)A(T). 



(17) 



According to Eqs. (|A15I) . (IA16[) and (|A17|) the coeffi- 
cients in these equations are given by 



t 



[3Mf - (M z \ H=0 y 



(18) 



fJ 



2 [a x + 2[i xz M*Y 2 [a y + 2f3 yz M^ 



(19) 



Here the angular brackets mean the averaging over the 
Fermi surface. JVq (k) is the angular dependent density of 
electronic states at the Fermi surface of the band t- The 
corresponding coefficients gjj and <?j| are obtained by the 
substitution f to i in Eqs. (fl~gj) and ([T^| . The function 
A(T) is 



A(T) 



2nTY — 

^ f,i... 



In 



n>0 



T' 



(20) 



where e is the energy cutoff for the pairing interaction. 
We assume here that it has the same value for both 
bands. 

The zero of the determinant of the system (TTTI) yields 
BCS like formula 



T = e exp 



1 

•911 



(21) 



where, the "constant" of interaction 




(22) 



is the function of temperature and magnetic field. Thus, 
the formula (1211) is in fact the equation for determina- 
tion of the critical temperature of the transition to the 
superconductig state T sc — T SC {H). 

Let us first consider the situation at H z = 0. Taking 
into account expressions (fT5|) and (IT^l) the equation (|2"T) 
can be rewritten as 



In- 



T 



1 

3|| 



-(T-T C ) 2 /(T), 



(23) 



where /(T) = const is slowly varying function of tem- 
perature. There are three different situations. 

(i) The Curie temperature is much larger than the cut- 
off energy, that is the band splitting is large in compar- 
ison with the width of region of effectiveness of pairing 
interaction. This case the solution of Eq. (l23l) obeys 
inequality 



(24) 



An increasing of the Curie temperature say as a func- 
tion of pressure shifts to the right the parabola in the 
right hand side of Eq. (|23|) . Hence, its intersection with 
logarithm in the left hand side is shifted to lower temper- 
ature. Thus, the critical temperature of transition to the 
superconducting state decreases with the Curie tempera- 
ture increasing as it is indeed the case in URhGe^ The 
opposite tendency takes place at decreasing of the Curie 
temperature accompanying by the increasing of the su- 
perconducting critical temperature as it is in UCoGe^^ 
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(ii) The outlined above behavior, however, has no 
proper description in the region where the Curie tem- 
perature decreases approching to the cutoff energy. This 
case the formal solution of Eq. (I2U1) is T sc = e = T c . It 
means that we are out applicability of BCS theory. So, 
the proper treatment of the problem in vicinity of inter- 
section of the T C (P) and T SC (P) curves is out of frame of 
applicability of developed approach. 

(iii) Finally one can consider the situation deeply in 
the paramagnetic region where the solution of Eq. (j2"3")l 
obeys inequality 



formulae 



T c < T sc < e. 



(25) 



Here, a displacement of parabola to the left shifts solution 
of Eq. (|23|) to the left. Hence, the Curie temperature 
decrease with pressure is accompanied by the decreasing 
of the temperature of superconducting transition. 

Let us turn now to the situation with H z ^ 0. In the 
ferromagnetic state the external field in the direction of 
the spontaneous magnetization drastically suppresses the 
superconducting state. Indeed, on the one hand the field 
dependence of the band density of states leads to the tiny 
shift in critical temperatue ST SC /T SC ~ ±^^lnj^M 
On the other hand, the increase of the magnetization un- 
der the influence of the parallel external field effectively 
suppress the coupling " constants" . For instance, at tem- 
peratures well below the Curie temperature the zero field 
magnetization is almost temperature independent. And 
the right hand side of formula (|2ip is temperature in- 
dependent as well. Then the critical temperature field 
dependence 



T SC (H Z 



eexp 



1 



9\\{H Z ) 



(26) 



is determined by the field dependence of "constants" 
of interaction at, ... which are expressed through the 
field dependent magnitude of magnetization. The mag- 
netic moment M z ° under magnetic field ~ 0.5 Tesla, 
which is of the order of the upper critical field at zero 
temperature, acquires ~ 1.3 times increase in respect to 
its zero field valued. According to Eq. (fT8|) it leads to 
the 4 times decreasing of the coupling constants and 
<jr* ! The decreasing of the constants and g^ (see Eq. 
(fl9)0 is not so impressive, but also certainly takes place. 
So, the magnetic field parallel to the spontaneous mag- 
netization effectively suppresses the pairing in the ferro- 
magnetic state. 

In the paramagnetic state in presence of the external 
field the coefficients g-\ ^ g± and g^ ^ gi and given by 



V z J z xx (k 2 x N 0t (k)) 
2[a z + 6f3 z M 2 } 2 ' 



v yJ y x (kiN 0t (k)) 



(27) 



(28) 



2 [a x + 2p xz M*Y 2 [ay + 2j3 yz M*Y 
where the magnetization has the field proportional value 
(see Eq. (|A18p ). Thus, above the Curie temperature the 
external magnetic field along the easiest magnetization 
axis suppresses the superconducting state as well as it 
does in ferromagnetic state for the field orientation par- 
allel to spontaneous magnetization. 



IV. MAGNETIC FIELD PERPENDICULAR TO 
MAGNETIZATION AXIS 

For the field orientated in a perpendicular to the spon- 
taneous magnetization direction say along b-axis it is 
natural to choose the spin quantization axis along the 
direction of the total magnetic field hz + H y y. Then the 
normal state matrix Green function is diagonal 



G n — 



G 1 ' 
G^ 



where 



G t,i - 



1 



iu n - £ k ± fi B J h 2 + Hi 



(29) 



(30) 



The potential of interaction found in the Appendix 
should be rewritten in the new coordinate frame 

7y(k,k') = R a Qrr^ x (k,k')5 im - £>L( k > k ') Rmj, 

(31) 



where 



1 
R = | cos tp — sin ip | , 

sin ip cos tp 



(32) 



is matrix of rotation around x direction on the angle given 
by tamp = H y /h, and R m j is the transposed matrix. 

The system of self-consistency equations acquires the 
following form 



At(k) - -T J2 n Ek' { [ D z ( k > k ') cos2 V + D y ( k > k ') sin2 f] G\GlAHk') 
+ [£>i(k,k')cos 2 (^ + (D^(k,k') - J D 2 L (k,k'))sinV] GtG^A i (k') 
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+i 



[D^(k, k') - D^(k, k')] sin ip cos tp (g\G^ + G|g|) A^(k')} , 



(33) 



A±(k) = -TJ2„ Ek< { [^-( k > k ') cos 2 + (l£(k, k') - ^(k, k')) sin 2 tp] G^G|At(k') 
+ [D^-(k,k')coB a vJ + I>iJ-(k,k / )Hin a v] G^A-^') 
+i [£>^-(k, k') - £>^(k, k')] sin cos tp {g\g\ +G{Gl) Ati(k')} , 

At4-(k) = -iTE„E k ' {~i [^(k,k') -^(k,k')] sin2p [GlG|At(k') + G^A^')" 
[(£>^(k, k') - D^(k, k')) cos 2 p + (Df (k, k') + D^(k, k')) sin 2 (g\G 1 2 + G\g\^ A^(k')} . 



(34) 



(35) 



Unlike to the case of parallel field here all the compo- 
nents of the order parameter are coupled each other. In 
principle one can find the solution of the whole system. 
However, in view of large band splitting one may neglect 



the last term containing the products GjG^+Gf G| in the 
equations ([33]) and ([34]) . Then the equations for thcA^ 
and A^- are decoupled from the equation (|3"5"j). The sys- 
tem acquires the same form 



At(k) = -TJ2 n J2 k , {[^(k,k')cos^ + ^(k,k')sin 2 ^] G^At(k') 
+ [Ui(k, k') cos 2 tp + (l£(k, k') - Z>^(k, k')) sin 2 tp] G{ G^(k')} , 



(36) 



A^(k) = -T£„ £ k , { [^-( k > k ') cos 2 <p + (D^(k, k') - D^K k')) sin 2 tp] G\g\/^ ik!) 
+ [L>^(k,k')cos 2 ^ + J D^(k,k')sin 2 ^] GfG^(k')} , 



(37) 



as it had in the case of parallel field given by Eqs. ([10]) 
and ([TT|) . Hence, assuming again that the highest critical 
temperature corresponds to solution given by Eqn. ([14")) 
and performing all necessary integrations and summation 
in Eqs. (l36l) and ([37)) we come to the system of algebraic 
equations 



n 



^(sV + sV)^) 



(38) 



According to Eqs. (1A"28)) . (|A"29)| . (fA30)) the coefficients 
in these equations are given by 



t _ VzlZx(*%No ( k )) 
' VWL(*^( k )> 



COS + 



^7j„(*«(k)> 



2 fa,. 



2 sin tp, 



(39) 



-t 



v 2[a, + 2/3 xz M|] 2 

/ V x l x xx (k 2 x N^k)) 
2[a x + 2p xz M%] 2 



V v ^ x (k 2 x N^))\ cQs2(fi 



2 [ay + 2p yz M%] 
32&M* 



sin 2 tp (40) 



r 



The corresponding coefficients and gjj are obtained by 
the substitution t to | in Eqs. and (|4H . 

The zero of the determinant of the system (fT71) yields 
BCS like formula 



T = e exp 



g± 



where, the "constant" of interaction 



9-L 




(41) 



(42) 



is the function of temperature and magnetic field. 

The easiest way to follow up the critical temperature 
field dependence is to consider temperature region well 
below the Curie temperature where the zero field mag- 
netization is almost temperature independent. And the 
right hand side of formula (|4ip is temperature indepen- 
dent as well. Then the critical temperature field depen- 
dence 



T sc (H y ) = eexp I - 



1 



g±(H y ) 



(43) 
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is determined by the field dependence of "constants" of 
interaction g\,g\, ... originating from the field dependent 
trigonometric factors 



cos 2 ip 



h 2 



H 2 '' 

v 



sin 2 ip 



HI 



H 2 

v 



(44) 



and the field dependent magnitude of magnetization 
M z (H y ) to be determined from Eqs ([MB and (|A"22j) . 

The trigonometric factors are changed with H y on the 
scale of exchange field. Hence, if the M z (H y ) is just 
slightly decreased at H y ~ H C 2 V <C h, as it is in URhGe in 
the region of moderate fields^, where M z (H y — 0) is big 
enough, then the critical temperature slowly decreases 
due to field dependence of trigonometric factors; 



T ac (H y ) - T sc (0) 
T sc (0) 



h 2 



(45) 



Obviously, this corresponds to small paramagnetic sup- 
pression of superconductivity have been considered in the 
paper—. 

Another situation takes place in UCoGe. Here the 
magnitude of zero field magnetization M z (H y = 0) is 
significantly smaller. Apparently, the field induced de- 
crease of M z (H y ) in this material is faster—, that leads 
to the effective increasing the "constants" of interaction 
g^, g^ despite of its decrease due to the trigonometric 
factors. This mechanism explains the effect of upward 
curvature in the temperature dependence of the upper 
critical fields in a and b directions observed in UCoGe^~— . 
At the same time the phenomenon of S shape tempera- 
ture dependence of the upper critical field in b-direction 
apparently related with significant increase of effective 
rnass^ observed in field interval from 5 to 15 Tesla. Out 
this interval for field in 6-direction and for any field in 
a-direction the average value of effective mass decreases 
with field increase and cannot be responsible for the ef- 
fect of the stimulation of superconductivity. 



V. CONCLUSION 

We have studied the magnetic field dependence of pair- 
ing interaction induced by magnetic fluctuations in the 
ferromagnetic superconducting compounds URhGe and 
UCoGe with orthorhombic crystal structure. 

For the field orientation along the spontaneous mag- 
netization there was demonstrated the effect of suppres- 
sion of fluctuations and, hence, of the superconductiv- 
ity by magnetic field. The leading role plays here the 
suppression of longitudinal fluctuations determining the 
constants g^ and g^ given by Eq. (fT5)l . that corresponds 
strongly anisotropic Izing like magnetism in uranium 
compounds^" 2 ^. 

For the field directed perpendicular to the spontaneous 
magnetization the field dependence of the constants of in- 
teraction is determined by interplay of paramagnetic sup- 



pression of superconductivity which is significantly weak- 
ened in the ferromagnetic superconductors^ and the ef- 
fect of stimulation of superconductivity by means the in- 
creasing of longitudinal magnetic fluctuations magnetic 
field given by Eq. ([39]). 
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Appendix A: Magnetic susceptibility of 
orthorhombic ferromagnet 

URhGe and UCoGe are the orthorhombic ferromag- 
nets with spontaneous magnetization oriented along c 
crystallography axis. At the temperatures below the 
Curie temperature and in the absence of magnetic field 
the c component of magnetization has a finite value. Here 
we derive the magnetic susceptibility of a ferromagnet 
with orthorhombic symmetry. It determines the pair- 
ing interactions due to spin fluctuations in such type of 
materials. To begin we write the Landau free energy of 
orthorhombic ferromagnet in magnetic field H(r) 



T = J dV(F M + F v ), 



(Al) 



where 



F M = a z M 2 + a y M 2 + a x M 2 
+p z M z A + p yz M 2 M y 2 + (3 XZ M Z 2 M X 2 - MH, (A2) 

and 

x dM x 8M x y dM y dM y z dM z 0M z 

i*v = 7 1" 7 + 7 

y dxi dxj v dxi dxj 13 dxi dxj 

(A3) 

is the density of gradient energy. Here the x,y,z are di- 
rections of the spin axes pinned to a, b, c crystallographic 
directions correspondingly, 



a z = a z0 (T - T c ), 



(A4) 



a x > 0, a y > and matrices 7?-, where p = x,y, z, have 
the form 



7L 

7 p y 



(A5) 
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Starting this point we shall separately consider the fer- 
romagnet susceptibility under stationary magnetic field 
parallel either to the magnetization direction or parrallel 
to one of the perpendicular crystallographic directions; 25 



1. Magnetic field parallel to magnetization axis 

So, let us take the magnetic field in the form 

H(r) = H z z + SH{r), (A6) 

where |(5H(r)| <C H z . By variation of the functional 
Eq. (|Al|) in respect to M z we come to 



where V z > is the constant of interaction. 

In similar manner one can obtain the expressions for 
the interaction functions in two other directions 



D|(k,k')« VxlfM ' ] 2 



and 



4( k > k ') 



2 [a y + 20y,M*\ 



(A15) 



(A16) 



2a z M z +4f3 z M*+2(3 yz M z M 2 -<y* 



d 2 M z 



where we have introduced the constants of interaction 
V x ,V y . Taking into account Eq. (|A10[) one can rewrite 
the expression for D z (k, k') as follows 



v '« dxidxj 



-H z -8H z {v)={) 
(A7) 



We shall search for the solution of this equation as 

M z (r) = M z + SM z (r), (A8) 

where M z is the solution of the space homogeneous prob- 
lem 



2a z M z + lfi z M° - H z =0. 



(A9) 



Here we omit the term 2j3 yz M z M 2 in view of zero coor- 
dinate independent value of magnetization 

My = 

in y direction. In absence of field the magnetization takes 
value 



{M z \ H=a) 2 = -g. 
For the 8M Z we have equation 



(A10) 



d 2 8M 

2a z SM z + 12(3 Z M 2 5M Z - 7 ?. = SH Z (r) (All) 

The solution of this equation in the k space is 

5M z (\l) = X z(k)8H z (k), (A12) 

where 



1 



2^ + 12^1+7?.^ 



(A13) 



The function determining the pairing interaction in 
triplet channel 19 is obtained from the subseptibility 
X z (k — k') as its part which is odd function of both of its 
arguments 



D\ (k, k') = -Yl [ Xz ( k - k') - Xz (k + k') 



2 [a z + 6/3 z M 2 } 



r, (A14) 



Dl(k,k') 



3 2 Z [3M 2 -(M Z \ H=0 ) 2 Y 



(A17) 



In the paramagnetic state the functions Z? z (k,k'), 
D x (k, k') and D y (k, k') are given by the same Eqs. 
(IA14P , (|A15I) , (|A16|) where the equilibrium value of mag- 
netization is proportional to the external field 



M, 



2a z 



(A18) 



Here, one must take in mind that the magnetization 
growing up at Curie temperature when a z (T) — >• is 
limited by the nonlinear term in Eq. (|A9|) . 



2. Magnetic field perpendicular to magnetization 

axis 

So, let us take the magnetic field in the form 

H(r) = H y y + 5H(r), (A19) 

where |<5H(r)| H y . Stationary vaue of magnetization 
is determined by the system of equations 

2a y M y + 2f3 yz M 2 z M y - H y = 0, (A20) 



a z + 2(3 Z M 2 + f3 yz M 2 = 0. 



(A21) 



The field induced stationary magnetization along b- 
direction is 



2(a y + /3 yz M 2 



(A22) 



Substituting this value back in the eqn. (|A2I) we obtain 
at f3 yz M 2 /a y < 1, that is certainly true not so far from 
the Curie temperature, 

F = a 0z (t~T c + P^I)m z 2 + /3 z M z 4 . (A23) 
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Hence, the Curie temperature 



(A24) 



is suppressed by the magnetic field oriented along b-axis. 
This type of behavior was observed in UCoGe£ The mag- 
netization along z-direction is also decreased 



Mi = 



a z0 (T c -T) PyzHy 



2/3 2 



(A25) 



Far from the Curie temperature, where this particular 
formula is not valid, the decrease of z-component of mag- 
netization determined by Eqs. (|A21[) and (|A22[) with 
growth H y still takes place. 

Following procedure of the previous section of Ap- 
pendix we obtain the susceptibility along z didection 



X*(k) = 



1 



2a z + + 20 VX M* + 7?.fcjfcj ' 

that making use Eq. (IA21[ ) can be rewritten as 

1 



X*( k ) = o - 



(A26) 



(A27) 



The corresponding function determining the pairing in- 
teraction is 



£^(k,k') 



32%M* 



(A28) 



By comparison with Eq. (|A17p in the absence of field we 
naturally come to = D\. 

The expressions for the interaction functions in two 
other directions are 



z^(k,k')« Vxlf ^ 2 

2 [a x + 2f3 xz M?} 2 



and 



^(k,k') 



2W, 



(A29) 



(A30) 



They have the same form as for the external field directed 
along the spontaneous magnetization. However, for the 
field orientation perpendicular to the spontaneous mag- 
netization M z decreases with field increase that leads to 
the growth the pairing interaction and increase the tem- 
perature of superconducting phase transition. 
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